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Abstract 


During  the  Shipboard  Integration  of  Sensor  and  Weapons  Systems  (SISWS)  Technology 
Demonstration  Project  (TDP),  an  algorithm  was  developed  to  automatically  generate  the  optimal 
chaff  solution  as  a  function  of  missile  initial  range  and  velocity,  wind  velocity,  and  engagement 
geometry.  The  chaff  solution  algorithm  was  developed  using  customized  score  functions  and 
simple  models.  The  SISWS  chaff  seduction  algorithm  was  improved  by  taking  a  more  physical 
approach  in  the  selection  of  chaff  solutions,  including  a  monopulse  antenna  model.  The 
document  starts  with  the  scenario  being  modeled,  followed  by  models  for  chaff,  intersect  ratio, 
monopulse  seeker,  and  probability  of  soft-kill. 


Resume 


Dans  le  cadre  du  Programme  de  demonstration  de  technologies  (PDT)  sur  T  integration  de 
capteurs  et  de  systemes  d’armes  embarques  (SISWS),  un  algorithme  a  ete  elabore  pour 
determiner  automatiquement  la  composition  optimale  des  paillettes  en  fonction  de  la  portee  et  de 
la  vitesse  initiales  d’un  missile,  du  vecteur  vent  et  de  la  geometrie  d’engagement.  L’algorithme  de 
determination  du  largage  des  paillettes  a  ete  elabore  a  Taide  de  fonctions  de  caracterisation 
adaptees  et  de  modeles  simples.  L’algorithme  de  seduction  par  paillettes  du  SISWS  a  ete  ameliore 
grace  a  une  approche  plus  concrete  pour  la  selection  de  la  composition  des  paillettes,  notamment 
grace  a  un  modele  d’antenne  monopulse.  Le  document  debute  par  la  modelisation  du  scenario, 
pour  ensuit  presenter  les  modeles  des  paillettes,  des  rapports  d’ intersection,  de  Tautodirecteur 
monopulse  et  de  la  probability  de  neutralisation  par  deroutement. 
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Executive  summary 


Improved  chaff  solution  algorithm 

Sylvain  Gauthier;  DRDC  Ottawa  TM  2008-324;  Defence  R&D  Canada  -  Ottawa, 
March  2009. 

During  the  Shipboard  Integration  of  Sensor  and  Weapons  Systems  (SISWS)  Technology 
Demonstration  Project  (TDP),  an  algorithm  was  developed  to  automatically  generate  the  optimal 
chaff  solution  as  a  function  of  missile  initial  range  and  velocity,  wind  velocity,  and  engagement 
geometry.  The  chaff  solution  algorithm  was  developed  using  customized  score  functions  and 
simple  models. 

The  SISWS  chaff  solution  algorithm  was  improved  as  described,  in  this  document,  by  adding  a 
monopulse  antenna  model  to  the  seeker,  and  by  using  a  more  physical  approach  to  the  selection  of 
solutions.  The  document  starts  with  the  scenario  being  modeled,  followed  by  models  for  chaff, 
intersect  ratio,  monopulse  seeker,  and  probability  of  soft-kill.  The  generated  chaff  solutions 
should  provide  more  accurate  responses  for  optimal  survivability. 

The  next  step  would  consist  of  analysing  the  efficiency  of  the  chaff  solutions  in  SADM  and  also 
comparing  the  results  with  the  SISWS  chaff  solution  algorithm.  After  that,  the  algorithm  can  be 
modified  further  to  include  the  tactic  of  combining  Continuous  Wave  (CW)  jamming  with  chaff 
solutions. 
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Sommaire 


Improved  chaff  solution  algorithm 

Sylvain  Gauthier;  DRDC  Ottawa  TM  2008-324;  Defence  R&D  Canada  -  Ottawa, 
Mars  2009. 

Dans  le  cadre  du  Programme  de  demonstration  de  technologies  (PDT)  sur  1’ integration  de 
capteurs  et  de  systemes  d’armes  embarques  (SISWS),  un  algorithme  a  ete  elabore  pour 
determiner  automatiquement  la  composition  optimale  des  paillettes  en  fonction  de  la  portee  et  de 
la  vitesse  initiales  d’un  missile,  du  vecteur  vent  et  de  la  geometrie  d’engagement.  L’algorithme  de 
determination  du  largage  des  paillettes  a  ete  elabore  a  l’aide  de  fonctions  de  caracterisation 
adaptees  et  de  modeles  simples. 

L’algorithme  de  composition  des  paillettes  du  SISWS  a  ete  ameliore,  tel  qu’indique  dans  le 
present  document,  grace  a  l’ajout  d’un  modele  d’antenne  monopulsee  dans  l’autodirecteur  et  a 
l’utilisation  d’une  approche  plus  concrete  pour  la  selection  de  la  composition  des  paillettes.  Le 
document  debute  par  la  modelisation  du  scenario,  pour  ensuit  presenter  les  modeles  des  paillettes, 
des  rapports  d’ intersection,  de  l’autodirecteur  monopulse  et  de  la  probability  de  neutralisation  par 
deroutement.  Les  compositions  de  paillettes  rendues  devraient  offrir  une  sensibilite  plus  precise 
en  vue  d’optimiser  la  capacite  de  survie. 

L’etape  suivante  consisterait  a  analyser  l’efficacite  des  compositions  de  paillettes  des  munitions 
de  destruction  atomique  speciale  et  de  comparer  les  resultats  a  l’algorithme  de  composition  des 
paillettes  du  SISWS.  Par  la  suite,  l’algorithme  pourra  etre  modifie  de  nouveau  pour  y  inclure  un 
brouillage  tactique,  ou  le  brouillage  a  ondes  entretenues  (CW)  precederait  le  largage  des 
paillettes. 
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1  Introduction 


During  the  Shipboard  Integration  of  Sensor  and  Weapons  Systems  (SISWS)  Technology 
Demonstration  Project  (TDP),  an  algorithm  was  developed  to  automatically  generate  the  optimal 
chaff  solution  as  a  function  of  missile  initial  range  and  velocity,  wind  velocity,  and  engagement 
geometry  [l]-[2].  An  initial  analysis  on  chaff  effectiveness  was  conducted  using  the  Ship  Air 
Defence  Missile  (SADM)  simulator  against  given  set  chaff  deploy  parameters  (launcher  number 
and  fuse  time).  Characteristics  of  the  most  successful  chaff  solutions  were  identified  from  the 
simulated  scenarios.  The  chaff  solution  algorithm  was  then  developed  using  customized  score 
functions  and  simple  models. 

Figure  1  shows  a  high  level  algorithm  developed  under  SISWS  for  the  generation  of  an  optimal 
seduction  chaff  solution.  The  parameters  of  the  scenarios  needed  to  generate  the  chaff  solutions 
are  first  extracted.  After  that,  a  set  of  candidate  solutions  for  chaff  launcher  and  fuse  time  is 
generated.  These  candidate  solutions  are  then  graded  against  customized  score  functions  for 
visibility,  miss  distance,  and  so  on.  The  best  solution  is  selected  by  finding  the  one  with  the 
highest  probability  of  capture  and  miss  distance. 

The  SISWS  chaff  seduction  algorithm  can  be  improved  by  directly  taking  into  account  the 
physics  of  the  problem  of  generating  the  chaff  solutions.  This  is  the  subject  of  this  document. 


Descript. 


Figure  1:  Overall  Seduction  Chaff  Solution  Algorithm 
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2  Modeling 


This  section  describes  the  models  used  in  the  improved  chaff  solution  algorithm.  The  main 
difference  in  modeling,  compared  to  the  SISWS  algorithm,  is  the  addition  of  a  monopulse 
antenna  model  for  the  missile  seeker.  A  more  physical  approach  was  also  used  in  the  selection  of 
chaff  solutions.  This  section  starts  with  the  scenario  being  modeled,  followed  by  models  for 
chaff,  intersect  ratio,  seeker  monopulse  antennas,  and  probability  of  soft-kill. 

2.1  Scenario 

An  antiship  missile  (ASM)  is  flying  at  constant  speed  and  constant  height  directly  toward  a  ship 
and  its  seeker  is  locked  onto  the  ship.  The  seeker  range-gate  width  (RGW)  is  centered  on  the 
ship.  The  ship  is  assumed  to  have  uniform  Radar  Cross  Section  (RCS)  of  37dB  for  all  look 
angles. 


2.2  Coordinate  system  transformations 

The  position  of  the  missile  in  the  ship  coordinate  system  can  be  described  in  Cartesian  or  polar 
coordinates  as  shown  in  Figure  3.  The  y  axis  is  along  the  ship  heading  and  the  x  axis  is 
perpendicular  to  the  y  axis  toward  the  right.  In  polar  coordinates,  the  position  of  the  missile  is 
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defined  by  its  range  to  the  ship  and  the  direction  angle  relative  to  the  y  axis.  The  coordinates  of  a 
point  p  can  be  transformed  from  one  system  to  the  other  as  follows. 


X  =  Rp  cos (6p );  y  =  Rp  sin(0p ) 


(1) 


0  =  tan  (xp/yp) 


(2) 


The  position  of  the  missile  can  also  be  described  into  the  RGW  coordinate  system  which  is 
centered  on  the  ship  but  with  its  axis  as  shown  in  Figure  4  and  Figure  5.  The  down  range  axis 
denoted  by  “d”  is  along  the  line  joining  the  missile  and  the  ship,  and  is  oriented  away  from  the 
missile.  The  cross  range  axis  denoted  by  “c”  is  perpendicular  to  the  d  axis  and  positive  to  the 
right.  The  coordinates  of  a  point  (xp,  yp)  can  be  transformed  into  the  RGW  coordinates  as  shown 
in  equation  (3). 

C  =  -  XP  CO<6 Missile  )  -  ^Missile  d  =  '  Missile  )  -  ^<6 Missile  )  (3) 


Missile 

✓ 


Figure  3:  Cross  range  and  down  range  coordinates  in  the  RGW  system 
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9  ®  Missile-  0° 


Figure  4:  Cross  range  and  down  range  coordinates  in  the  RGW system  (6missiie  =  0°) 


c 


d 


(Xp,  yp) 
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^Missile- 


Figure  5:  Cross  range  and  down  range  coordinates  in  the  RGW  system  (6missiie  =  90°) 


4 


DRDC  Ottawa  TM  2008-324 


2.3  Chaff  trajectory 

2.3.1  Range  height  profile  of  chaff  round 


The  range  and  height  of  the  chaff  round  is  a  function  of  fuse  time  (Tf)  as  shown  in  Figure  6  (for 
three  elevation  barrels)  [1].  Ignoring  the  rolls  and  pitches  of  the  ship,  the  chaff  barrel  elevation  is 
always  30°.  Table  1  shows  the  extracted  values  of  chaff  range  and  height  for  the  30°  barrel 
elevation. 

The  interesting  point  here  is  the  fact  that  chaff  rounds,  if  not  fused  fall  back  to  launcher  height 
after  about  19.5  seconds.  This  limits  the  range  of  possible  fuse  time  solutions  from  zero  to  19.5 
seconds.  Indeed,  chaff  launcher  height  is  about  5  meters  above  sea  level.  The  fall  rate  of  chaff  is 
typically  between  0.6to0.9m/s  [1].  Hence  a  chaff  cloud  would  fall  to  the  sea  surface  about  6 
seconds  before  it  blooms  to  90%  of  its  full  value. 


Figure  6:  Chaff  range  and  height  as  a  function  of fuse  time  Tfand  barrel  elevation 
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Table  1 :  The  range  and  height  of  chaff  round  as  a  function  of  fuse  time  (30°  in  elevation) 


TOF  (sec) 

Range  (m) 

Height  (m) 

TOF  (sec) 

Range  (m) 

Height  (m) 

0 

0.0 

0 

10 

1,659.1 

420 

1 

90.9 

50 

11 

1,795.5 

410 

2 

295.5 

120 

12 

1,931.8 

390 

3 

500.0 

210 

13 

2,045.5 

360 

4 

681.8 

270 

14 

2,181.8 

320 

5 

863.6 

320 

15 

2,318.2 

280 

6 

1,045.5 

360 

16 

2,431.8 

220 

7 

1,204.5 

390 

17 

2,522.7 

160 

8 

1,363.6 

410 

18 

2,659.1 

110 

9 

1,522.7 

420 

19 

2,750.0 

60 

19.5 

2,795.5 

0 

A  mathematical  model  of  the  chaff  round  is  range/height  profile  in  Table  1  was  produced  by 
using  a  Matlab  function  called  polyfit,  with  the  degree  set  to  six.  The  resulting  polynomial 
coefficients  that  fit  these  data  define  the  chaff  range  /  height  profile,  as  shown  below.  The 
range/height  profiles  of  the  chaff  round  obtained  from  equations  (4)  and  (5)  compare  very  well 
with  the  original  data  as  shown  in  Figure  7. 


rchaff  =  0.003t6  -0.0225t5  +0.589U4  -7.4036?3  +41.8246?2  +91. 174k  (4) 

Khaff  =  0.000k6  -0.0077f5  +  0.204k4  -2.6225t3  +11.3826/2  +52.4026f  (5) 


Once  the  chaff  round  explodes,  its  dimension  reaches  100m  along  the  launch  axis  [1].  The 
exploding  rounds  will  be  spread,  centered  on  its  fuse  position.  Hence,  the  chaff  round  should  not 
be  fused  before  traveling  at  least  50  meters.  Using  equations  (4)  and  (5),  this  minimum  distance 
corresponds  to  a  fuse  time  of  0.41  seconds.  Typical  tactics  have  a  minimum  fuse  time  of  0.75s, 
which  corresponds  to  a  round  distance  of  99.5  meters. 
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Chaff  range  and  height  for  fuse  time  Tf 
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Figure  7:  Chaff  range  and  height  vs  fuse  time  Tfand  barrel  elevation 


2.3.2  Round  position  at  fuse  time  relative  to  ship  x,  y  coordinates 

The  positions  of  each  chaff  launcher  are  given  in  Table  2  and  shown  in  Figure  8  [1].  The  initial 
position  of  the  chaff  round  at  fuse  time  in  x,  y  coordinates  are  then  given  in  equations  (6),  (7)  and 
(8). 


x chaff  init  (L# ,  Tf )  —  xL#  H  rchaff  (Tf )  *  cos(90  6az )  (6) 

y  chaff  init  (L#-,Tf)  =  yL#  H  rchaff  (Tf)*  sin(90  -  6az )  (7) 

Z  chaff  init  (L#  ,  Tf  )  =  ZL#  H  hchaff  (Tf  )  (8) 
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Table  2:  Position  of  each  launch,  in  shi] 

3  coordinates 

Launcher  # 

x(m) 

y  (m) 

z(m) 

0az  (degree) 

1 

8 

15 

5 

30 

2 

-8 

15 

5 

-30 

3 

8 

-15 

5 

110 

4 

-8 

-15 

5 

-110 

Height  of  launcher  is  assumed  to  be  5  meters  above  sea  level. 


2.3.3  Cloud  Drift 

We  assume  that  the  chaff  cloud  at  fuse  time  is  centered  on  the  round  position  at  Tf  (mass  center). 
Each  part  of  the  cloud  then  drifts  along  the  wind  direction  at  the  same  speed.  Annex  A  shows 
distance  traveled  by  a  chaff  cloud  in  function  of  wind  speed  and  elapsed  time.  Chaff  cloud 
typically  has  a  fall  rate  between  0.6  to  0.9  m/s.  Let  us  define  dt  as  the  elapsed  time  after  fuse 
time.  Then  the  position  of  the  chaff  center  is  given  by  equations  (9),  (10)  and  (11). 


x 


chaff  dt y 


(L# ,  Tf )  —  xchaff  init  ( L # ,  Tf )  H  Vwind  x*  dt 


(9) 


y  chaff  dt  (L#  ?  Tf )  —  ychaffinit  (L# ,  Tf )  H  Vwind  y*  dt 


(10) 
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^  chaff  dt 


(L#,Tf)  =  zchaffinit(L#,Tf) 


V  *  dt 

V  chaff  fall  Ul 


(11) 


Using  these  equations,  the  trajectory  of  the  chaff  cloud  for  all  acceptable  fuse  time  solutions  (0  to 
19.5  seconds)  and  launcher  number  can  be  plotted. 

The  cross  range  and  down  range  position  of  the  chaff  cloud  in  the  RGW  coordinate  system  can  be 
calculated  using  previous  equations  of  transformation  coordinates  resulting  in  equations  (12)  and 
(13). 


c  chaff  dt  (L# ,  Tf )  —  xchaffdt  cos  (6  Missile )  H  y  chaff  dt  (L#,Tf)  sin  (6Missile )  (12) 


d  chaff  dt(L#,Tf)  —  ychaffdt  (L#,Tf)  cos(6  Missile )  x  chaff  dt(L#,Tf)  sin  (0Missile)  (13) 


2.4  Intersect  Ratio 

This  subsection  calculates  the  amount  of  chaff  cloud  that  intersects  with  the  seeker  RGW  which  is 
centered  on  the  targeted  ship.  Figure  9  shows  an  example  of  chaff  cloud  center  trajectories  from 
launchers  1  and  4,  which  intersect  with  the  missile  seeker  RGW.  This  subsection  starts  by 
modeling  the  physical  extent  of  the  chaff  cloud  and  associated  trajectory  in  the  RGW  coordinate 
system.  The  amount  of  intersection  with  the  RGW  is  then  calculated  using  equation  (13). 
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Figure  9:  Ship,  missile,  wind  engagement  geometry. 


2.4.1  Cloud  extent  in  down  range 

The  chaff  cloud  is  first  centered  on  the  round  position  at  fuse  time  and  then  every  part  of  the 
cloud  is  drifted  by  the  wind  in  the  same  direction  and  speed.  The  extent  of  the  chaff  cloud 
relative  to  the  chaff  center  is  given  by  equations  (14),  (15)  and  (16). 


x 


chaff  _  extent 


L 


chaff  * 

2 


Wr 


sinC^-  Launcher  )  +  -|SL.  cos(0„ 

Launcher  ) 


(14) 


y  chaff  _  e 


Jdiaff_  -| ■ 

2 


W, 


cos  (6az 

Launcher  )  +  -^L*sin(0„ 

Launcher  ) 


(15) 


Z 


chaff  _  extent 


^ chaff 

2 


(16) 
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The  cross  range  and  down  range  extent  of  the  chaff  cloud  in  the  RGW  coordinate  system  is 
obtained  through  the  equations  of  coordinates  transformation,  (17)  and  (18). 


"  chaff  _  extent  ^  chaff  _  extent  C0S(^ Missile)  ^  y chaff  _  extent  Sin(^ Missile  ) 


(17) 


d  chaff  _extent=~  T  chaff  _extent  C0S(^  Missile)  ~  XCkaff  _6Xtent 


(18) 


The  down  range  values  of  the  chaff  cloud  extremities,  dchaff+  and  dchaff-,  in  the  RGW  coordinate 
system  are  given  by  equations  (19)  and  (20). 


d chaff  dt+  d chaff  dt  ^  d chaff  _  extent 


(19) 


d chaff  dt-  d  chaff  dt  d  chaff  _  extent 


(20) 


2.4.2  Intersection  with  the  RGW 


There  are  four  cases  where  the  cloud  of  chaff  can  intersect  the  RGW  as  shown  in  Figure  10. 
Using  equations  (19)  and  (20),  we  can  determine  if  one  of  these  intersects  cases  occurs.  The 
amount  of  chaff  in  down  range  that  intersects  with  the  RGW  can  then  calculated  using  the 
equation  associated  to  that  intersects  case,  (21),  (22),  (23)  or  (24). 


Case  1:  d( 


chaff  dt+  ' 


>  RGW/2  and  dchaff*-  is  inside  RGW 


[  chaff  dt 


(L#,T , , dtcase , )  — 


RGW 


-  d 


chaff  dt- 


(21) 


Case  2:  dchaff  dt-  <  -  RGW/2  and  dChaff  dt+  is  inside  RGW 


IchaffdtiLtt’Tf’dt  easel)  ~  d 


f  ?  U'L easel; 


1  chaff  dt+ 


+ 


RGW 


(22) 


Case  3:  dchaff  dt-  and  dchaffdt+  are  inside  RGW 

I  chaff  dt(Ltt>Tf>  dtCase  3)  =  d 


chaff  dt+ 


1  chaff  dt- 


(23) 
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Case  4:  dchaff  dt-  <  -  RGW/2  and  dchaff  dt+  >  RGW/2 

Ichaffdt(L#,Tf,dtcase  4)  =  RGW  (24) 


Figure  10:  Chaff  cloud  interseciont  cases  with  RGW 


The  ratio  (or  fraction)  of  chaff  intersects  the  RGW  is  given  in  equations  (25)  and  (26). 


ratio  dt 


( L#,Tf,dt )  = 


1  chaff  dt 


(. L#,Tf,dt ) 


1  chaff  _  extent 


(25) 


I  ratio  dBdt  (L# ,  Tj ,  dt)  —  1 0  log10  {I  ratio  dBdt  (L# ,  Tj- ,  dt))  (26) 

Using  these  equations,  the  fuse  time  solutions  that  result  in  the  chaff  trajectory  intersecting  the 
RGW  can  be  extracted  and  plotted. 
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2.5  Minimum  miss  distance 


The  minimum  value  of  miss  distance  in  cross  range  is  a  function  of  the  physical  extent  of  the  ship 
and  a  minimum  safe  distance  threshold  set  by  the  user.  Figure  1 1  shows  the  physical  extent  of  the 
ship  relative  to  the  RGW  coordinate  system.  It  is  found  that  the  minimum  value  of  miss 
distance  is  given  by  equations  (27)  through  (29). 


miss  mm 


=  max  (cn 


min  2  > 


(27) 


where 


"  miss  min  1 


W. 


ship 


:cos  {Gmit 


missile 


j  +  ^ship  * 


Sin(0,nissile) 


+  Thres 


mm  miss 


(28) 


^  miss  min  2 


W  L 

^L‘cos(0„„„)  +  -St*sin(0, 


missile  > 


+  Thres  min 


(29) 


and  Thres min  miss  is  the  minimum  cross  range  distance  between  the  missile  and  the  ship  for  the 

ship  to  be  safe.  Solutions,  where  the  chaff  miss  distance  in  cross  range  is  smaller  than  the  user 
minimum  value  of  miss  distance,  are  dismissed. 
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Figure  11:  Ship  dimension  in  x,  y  coordinates 

2.6  Chaff  RCS  Models 

The  model  for  the  chaff  RCS  as  a  function  of  elapsed  time  since  fuse  time  is  taken  from  reference 
[1].  In  this  model,  the  chaff  RCS  as  a  function  of  time  is  given  by  equation  (30). 


RCS 'chaff  dt  ~  RCSchaff  max  * 


l-e- 


(30) 


where  Kbloom  is  a  blooming  constant  which  is  set  as  2.2,  and  RCSchaffdt  is  in  dB.  Figure  12  shows 

the  value  of  RCS  as  a  function  of  elapsed  time  after  fusing.  The  chaff  takes  about  10  seconds  to 
reach  its  maximum  value,  which  is  relatively  long  considering  that  a  missile  is  coming  toward  the 
ship.  The  value  of  the  chaff  RCS,  relative  to  its  maximum  value,  is  given  by  equation  (31). 


RCS 


chaff  dtrel 


=  RCS 


chaff  dt 


-  RCS 


chaff  max 


(31) 
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Chaff  RCS  as  a  function  of  blooming  time 


Figure  12:  RCS  of  blooming  chaff  as  a  function  of  time 


2.7  Seeker  Antenna  model 

An  ASM  seeker  typically  uses  a  monopulse  antenna  system  to  track  the  targeted  ship.  Monopulse 
radar  consists  of  two  antenna  beams  offset  by  the  3dB  antenna  beamwidth  [3].  The  modeling  of 
the  whole  tracking  process  would  be  relatively  complex  and  would  require  a  lot  of  processing 
time.  A  simple  approach  is  used  in  this  document  to  model  the  probability  of  the  seeker  detecting 
the  chaff  and  being  captured  by  it.  This  subsection  starts  by  modeling  an  antenna  beam  pattern 
followed  by  the  monopulse  antenna  system. 


2.7.1  Antenna  beam  pattern 

In  the  far-field  the  electric-field  pattern  of  an  aperture  antenna  is  the  Fourier  transform  of  the 
electric  field  illuminating  the  aperture  [4]  to  [6].  For  uniformly  illuminated  one-dimensional 
aperture,  the  antenna  pattern  is  proportional  to  the  square  of  a  scaled  sine  function  [4].  The  sine 
function  is  given  by  equation  (32)  and  is  plotted  in  Figure  13. 


sinc(#)=  1,  (f=  0 

sin(^) 

— (p*  0 

K(p 


(32) 
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where  cp  is  in  radian.  The  3dB  beamwidth  (half  value  on  linear  scale)  of  the  sine  function  is 
equal  to  25.379  degrees.  The  sine  function  can  be  scaled  to  other  beamwidths  by  multiplying  the 
cp  argument  in  radians  with  the  ratio  of  the  3dB  beamwidth  of  the  sine  function  on  the  new 
beamwidth  as  shown  in  equation  (33). 


V scaled 


=  (p* 


thdBsinc 

VldBb  earn 


(33) 


The  normalized  antenna  gain  at  angle  cp,  relative  to  boresight,  is  given  by  equation  (34). 


G,el  = 


( 


smc 


r 


^  PsdBs 
VldBbec 


-|2 


/j 


(34) 


where  (f3dB^mc  is  the  3dB  beamwidth  of  the  non-scaled  sine  function  and  is  equal  to  25.379  in 

degrees.  (f3dBheam  is  the  3dB  beamwidth  of  the  antenna  beam  to  be  modeled.  Figure  14  shows  the 

antenna  pattern  using  the  square  of  the  scaled  sine  function  for  a  3dB  beamwidth  of  5°.  Annex  A 
shows  values  of  the  arc  sustained  by  a  beam  of  4  degrees  in  function  of  range.  The  length  of  the 
arc  sustains  by  the  seeker  main  beam  having  a  4  degrees  beam  will  be  about  700m. 
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Square  of  sine  Amplitude  of  sine 


Plot  of  sine  function 


Square  of  the  sine  function  in  function  of  theta  in  degree 


Figure  13:  Plot  of  the  sine  function  and  its  squaret 
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ANtenna  gain  pattern  in  decibel  for  a  3dB  beamwidth  of  5  degrees 


Figure  14:  An  example  of  antenna  gain  pattern  in  decibels  for  a  3dB  beamwidth  of  5  degrees 

2.7.2  Monopulse  antenna  system 

A  monopulse  antenna  system  uses  two  identical  antenna  beams  offset  by  an  angle  equal  to  the 
3dB  of  the  antenna  beam  as  shown  in  Figure  15.  From  the  figure,  we  have  equations  (35)  and 
(36) . 

0A=0»-?y  (35) 

(36) 


where  0A  is  the  boresight  direction  of  antenna  A;  0B  is  the  boresight  direction  of  antenna  B;  0M  is 
the  boresight  direction  of  the  monopulse  antenna  system;  and  03dB  is  the  3dB  antenna  beamwidth. 
The  chaff  angle  relative  to  the  boresight  of  antenna  A  and  B  are  given  by  equations  (37)  and  (38). 
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(37) 


^  chaff  _  A  ~  &  chaff  @  A 


@ chaff  _B  ~  &  chaff  ^  B 


(38) 


The  antenna  gain  at  the  chaff  position  for  each  antenna,  relative  to  their  boresight,  is  given  by 
equations  (39)  and  (40). 


chaff  _ A 


smcf  6  *  -  &3dBsinc 

SHIM  Uchaff  _A 


12 


V 


6 


3dBbeam 


(39) 


'chaff  _B 


f 


sine 


v 


p  *  ^  ^3dB sine 

° chaff  _B  q 


~|2 


3  dBbeam 


(40) 


where  Gchaff_A  and  Gchaff_B  are  the  antenna  gain  for  the  chaff  located  at  0chaff  from  antenna  A  and  B 
relative  to  their  boresight  directions. 
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Figure  15:  Monopulse  antenna  system 


2.7.3  Probability  of  detection 

The  probability  of  detecting  the  chaff  is  a  function  of  the  chaff  position  relative  to  the  monopulse 
boresight.  The  probability  of  detecting  a  target  for  each  antenna  channel  can  be  calculated  using 
the  equation  of  Neuvy  [7],  (41). 


^=10- 


where  u  — 


a  log 


/ln(2)x 

p» 


I  Hfi 


n2J3  SNR 


(41) 
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The  constant  a  and  [3  are  functions  of  the  Swelling  type  of  target  detected.  For  a  chaff  target,  (3 
would  be  equal  to  1.  In  our  case,  the  Signal  to  Noise  Ratio  (SNR)  cannot  be  calculated  since  the 
seeker  characteristics  are  not  known.  However,  the  antenna  seeker  is  assumed  to  be  locked  onto 
the  center  of  the  ship.  Hence,  the  probability  of  detection  should  be  equal  to  one  when  reflected 
power  is  greater  or  equal  to  the  power  received  from  the  ship.  As  the  ratio  of  chaff  to  ship  power 
decreases  below  unity,  the  probability  of  detecting  the  chaff  will  decrease  until  it  approaches 
zero.  Using  the  same  form  of  equation  as  Neuvy,  the  probability  of  detecting  the  chaff  can  be 
expressed  as  in  equation  (42). 


d  chaff 


10 

1 


P chaff  /  ship 


•  p  <  o 

I1  chaff  /ship  ^  u 

•  P  >0 

V  chaff  /ship  ~  u 


(42) 


where  Pchaff/ship  is  the  maximum  ratio  of  received  power  from  the  chaff  and  the  ship.  The  1  in 

the  exponent  was  added  to  get  a  probability  of  detection  of  unity  when  the  power  ratio  is  1 . 
Figure  16  shows  the  modeled  probability  of  chaff  detection  as  a  function  of  the  received  power 
ratio,  in  dB,  between  the  chaff  and  the  ship.  The  relative  antenna  gain  at  the  ship  location  is  3dB 
lower  than  on  the  boresight  of  each  antenna.  Hence,  the  probability  of  chaff  detection  decreased 
rapidly  beyond  the  3dB  beamwidth  of  both  antennas.  This  is  a  realistic  feature  of  a  monopulse 
system. 
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Probability  of  detecting  the  chaff  as  a  function  of  the  ratio  of  received  power  from  chaff  and  ship,  in  dB 


Figure  16:  Probability  of  detecting  the  chaff  when  the  monopulse  system  is  locked  onto  the  ship 


The  power  reflected  back  from  the  ship  and  the  chaff  into  each  antenna  channel  is  proportional  to 
the  RCS,  propagation  factor,  antenna  gain  relative  to  boresight  and  intersect  ratio,  as  shown  in 
equations  (43)  and  (44). 


(43) 


PshiP(?)  ®  c  ship (0  FshiP (0  Grel  ship(t)  I ship(t)  (44) 

where  C  is  RCS,  F  is  the  propagation  factor,  Grel  is  the  antenna  gain  at  chaff  position  relative 
to  boresight,  and  I  is  the  intersect  ratio  with  the  missile  RGW.  The  ratio  between  the  received 
power  from  the  chaff  and  the  ship  is  then  given  by  equation  (45). 
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(45) 


1  chaff  I  ship 


»«)  =  ■ 


O chaff  (0  F chaff  (0  ^ 


chaff v 


7  re/  chaff  max 


(0  I  chaff  (0 


^ ship  (0  ^ shiv  ^ rel  shiv  (0  ^ shiv  CO 


ship 


rel  _  ship ' 


1  ship " 


The  ship  is  located  on  the  monopulse  antenna  system  boresight  and  centered  on  the  RGW. 

Hence,  the  relative  antenna  gain  is  0.5  and  the  intersect  ratio  is  1.0.  Ignoring  propagation  effects 
such  as  multipath  we  have  equation  (46). 


‘  chaff  /ship 


*(0  =  - 


chaff 


it)  G 


rel  chaff  max 


(0  I  chaff  (0 


0-5 


(46) 


This  generation  of  chaff  solutions  does  not  consider  missile  height.  This  is  not  too  important 
since  the  seeker  is  locked  onto  the  ship,  i.e.,  it  is  within  the  seeker  beamwidth.  The  missile 
height  will  become  important  when  the  propagation  factor  is  included  in  the  model. 

2.7,4  Probability  of  capture 

An  analysis  of  the  monopulse  response  (A/£  channel),  as  the  chaff  moves  away  from  the  ship,  was 
done  using  the  output  of  the  monopulse  antenna  channels  (see  Annex  B).  It  was  found  that  the 
equivalent  RCS  center  will  be  tracked  until  the  monopulse  system  switches  to  one  of  the  targets. 
The  monopulse  tracker  will  come  back  to  the  ship  whenever  its  RCS  is  greater  or  equal  to  the 
chaff  RCS.  Also,  the  probability  of  recapture  by  the  ship  decreases  as  the  chaff  RCS  exceeds 
the  RCS  of  the  ship.  Based  on  previous  findings,  a  simple  model  for  the  probability  that  the  chaff 
will  stay  captured  until  the  end  is  given  by  equation  (47). 


p  =l-e 

1  capt  1  e 


and  ° 'chaff  ship  dB  >  0 


(47) 


When  the  chaff  RCS  exceeds  the  ship  RCS  by  3dB,  the  probability  of  chaff  staying  captured  until 
the  end  is  close  to  95%,  which  seems  realistic  [1].  Figure  17  shows  the  probability  of  chaff 
capture  as  a  function  of  reflected  power  relative  difference  based  on  equation  (47). 
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Probability  of  chaff  staying  captured  until  the  end  of  TTG 


Figure  1 7:  Probability  of  chaff  capture  in  function  of  reflected  power  relative  difference 


For  each  intersection  point  within  the  RGW,  the  probability  that  the  chaff  will  be  detected  and 
stay  captured  is  given  by  equation  (48). 


Pdcapt(t)  =  Pd(t)*Pcapt 


(48) 


The  total  probability  of  chaff  being  detected  and  captured  within  the  RGW  is  given  by  equation 
(49). 

Ptot_d_capt  =  1 ' -  f  ,  0  “  Pd_capt  (0)  (49) 
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2.8  Choice  of  best  solution 


Chaff  can  be  launched  immediately  at  initial  Time  To  Go  (TTG)  or  can  be  delayed  by  steps. 
Chaff  solutions  having  the  larger  number  of  time  delays  before  launch,  provide  more  flexibility  to 
the  operator.  The  best  solution  is  then  chosen  among  the  group  of  solutions  (Tf,  L#)  having  the 
greater  number  of  possible  time  delays.  The  best  solution  from  this  group  is  the  one  that  has 
greatest  weight  for  miss  distance,  detection  and  capture  defined  in  equations  (50)  and  (51): 


wtot  =  0.75  wd  capt  H  0.25  w 


Rmiss 


(50) 


^ d  _ capt  _d  _ capt  ’ 


^ Rmiss 


1 


Kiss 

Min  _  Rmiss 


(51) 


where  wd  capt  is  the  weight  associated  with  the  total  probability  of  being  detected  and  staying 

captured,  and  wRmiss  is  the  weight  associated  with  the  miss  distance.  A  greater  weight  was  given 

to  detection  and  capture  since  it  must  occur  to  be  successful.  The  weight  for  miss  distance  is  also 
lower  since  only  the  solutions  with  miss  distance  above  safe  thresholds  are  kept. 

The  weight  associated  with  the  miss  distance  was  defined  assuming  that  beyond  a  certain  distance 
in  miss  distance  there  is  not  much  difference  in  the  outcome.  Figure  18  shows  an  example  of 
weight  associated  with  the  miss  distance.  The  minimum  miss  distance  was  set  to  100m  in  this 
case.  As  it  is  seen,  the  weight  associated  with  the  miss  distance  decreases  rapidly  when  the  miss 
distance  exceeds  the  minimum  value  and  then  approaches  unity  around  300  meters. 
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Weight  for  Rmiss 


Figure  18:  Weight  associated  with  miss  distance 
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3  Improved  algorithm 


3.1  Algorithm 

Table  3  shows  the  improved  chaff  solution  algorithm,  at  a  high  level.  It  can  be  divided  into  four 
parts:  possible  solutions,  geometric  criteria,  seeker  criteria  ,and  selection  of  best  solution. 

This  algorithm  starts  by  generating  trajectories  of  all  possible  chaff  solutions  including  a  set  of 
discrete  time  delays  before  launch  of  chaff  round.  The  chaff  round,  as  seen  in  section  2,  will  fall 
back  at  launcher  level  after  about  19.5  seconds.  This  maximum  value  of  fuse  time  for  chaff 
solutions  is  then  set  to  19.5s  (or  initial  Time  To  Go  (TTG)  if  smaller  than  that).  The  minimum 
value  of  Tf  was  set  to  0.42s  to  avoid  exploding  the  chaff  round  over  the  ship.  Indeed,  the  chaff 
round  travels  50  meters  during  that  time,  which  corresponds  to  half  the  length  of  the  exploded 
chaff  cloud.  In  the  past,  only  small  subsets  of  candidate  solutions  for  fuse  time  were  considered. 

After  that,  geometric  considerations  are  used  to  filter  out  chaff  solutions  where  chaff  clouds  do 
not  intersect  with  the  seeker  range  gate  and  where  miss  distance  is  below  a  given  threshold.  The 
range  gate  is  here  assumed  to  be  infinite  in  length,  perpendicular  to  the  missile  trajectory. 
Equations  developed  in  section  2  are  used  to  filter  out  unwanted  solutions. 

After  that  seeker  considerations  are  used  to  filter  out  chaff  solutions,  which  the  total  probability 
of  detection  and  capture  is  below  a  minimum  threshold.  This  step  takes  into  account  the 
difference  in  RCS  and  received  power  between  the  chaff  and  the  ship.  The  corresponding 
equations  developed  in  section  2  are  used  to  filter  out  unwanted  solutions. 

The  next  step  consists  of  selecting  the  best  solutions  from  among  these  final  solutions.  The  group 
of  solutions  (Tf,  L#)  having  the  greater  number  of  possible  time  delays  before  launching  is  first 
selected.  The  best  solution  of  this  group  is  then  selected  based  on  a  weighted  sum  for  total 
probability  of  detection  and  capture,  and  also  for  the  miss  distance. 
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Table  3:  High-level  chaff  solution  algorithm 


Possible  solutions 

1 .  Calculate  trajectories  of  all  possible  chaff  solutions 


Geometric  criteria 

2.  Keep  solutions  that  intersect  with  RGW 

3.  Keep  solutions  that  exceed  minimum  miss  distance 


Seeker  criteria 

4.  Calculate  relative  difference  in  reflected  power 

5.  Keep  solutions  that  exceed  minimum  probability  of  capture 


Select  optimal  solution 

6.  Select  solution  that  optimize  probability  of  soft-kill 


3.2  Results  example 

The  scenario  considered  in  this  subsection  includes  a  missile  flying  at  300m/s  coming  from  the 
head  of  the  ship  at  an  initial  range  of  12,000m.  The  antenna  beams  of  the  seeker  have  a  3dB 
beamwidth  of  4°  and  the  seeker  RGW  is  equal  to  100m.  The  wind  speed  is  20  knots  away  from 
ship  toward,  190°  clockwise  from  ship  heading.  The  minimum  miss  distance  is  set  to  50m.  The 
minimum  probability  for  detection  and  capture  is  set  to  60%.  The  maximum  RCS  of  chaff  is  set 
to  40dB  and  the  average  ship  RCS  is  set  to  37dB.  Figure  19  shows  the  trajectories  of  all  possible 
chaff  solutions  for  this  scenario.  The  grey  circle  indicates  the  position  of  the  chaff  round  at  fuse 
time.  The  chaff  then  drifts  along  the  wind  direction  until  the  total  elapsed  time  (Tf  plus  drift 
time)  is  equal  to  the  initial  missile  TTG.  The  chaff  trajectories  are  not  symmetrical  on  each  side 
of  the  ship.  This  is  due  to  the  relative  wind  direction  heading  of  190°.  Figure  20  and  Figure  21 
show  the  solutions  from  previous  figures  where  the  chaff  clouds  intersect  with  the  missile  range 
gate  and  the  miss  distance  above  a  minimum  value.  Chaff  center  positions  within  the  minimum 
miss  distance  are  shown  in  red  to  enhance  awareness.  As  can  be  seen,  the  number  of  solutions 
has  been  significantly  reduced.  Figure  22  shows  the  chaff  solutions  where  the  total  probability  of 
chaff  detection  and  capture  is  above  the  minimum  value  set  by  the  user.  Once  again  the  number 
of  chaff  solutions  was  reduced.  All  solutions  from  launchers  3  and  4  were  discarded.  This  makes 
sense  since  their  trajectories  intersect  very  little  with  the  RGW.  The  best  solution  is  shown  in 
green  and  the  probability  of  soft-kill  is  one  in  this  case.  The  black  dots  on  the  chaff  trajectories 
are  the  final  positions  when  chaff  launches  are  delayed  until  the  missile  reaches  initial  Range  To 
Go  (RTG)  minus  0,  1.5,3,  4.5  and  6km.  The  more  delays  you  can  have,  the  more  robust  the 
solution.  For  example,  the  best  solution  is  still  successful  even  if  the  chaff  is  launched  when  the 
missile  RTG  is  only  6km  rather  than  12km. 
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y  axis  (meter)  y  axis  (meter) 


Trajectory  of  chaff  for  all  fuse  time  as  Tf,  in  ship  coordinates. 


Figure  19:  Trajectories  of  all  possible  solutions  for  a  given  scenario 
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Figure  20:  Trajectories  of  all  solutions  that  intersect  RGW for  a  given  scenario 
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Figure  21:  Close  up  of  the  trajectories  of  all  solutions  that  intersect  RGW for  a  given  scenario 
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Figure  22:  Trajectories  of  all  solutions  that  exceed  minimum  probability  of  soft-kill 
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Another  scenario  considered  here  is  a  missile  flying  at  300m/s  coming  from  the  head  of  the  ship 
and  initially  at  a  range  of  4,700m.  The  antenna  beams  of  the  seeker  have  a  3dB  beamwidth  of  4° 
and  the  seeker  RGW  is  equal  to  100m.  The  wind  speed  is  20  knots  and  is  varied  from  0  to  360 
degrees  going  away  from  ship  clockwise  from  ship’s  heading.  The  minimum  miss  distance  was 
set  to  50m.  The  minimum  probability  for  detection  and  capture  is  set  to  60%.  The  maximum 
RCS  of  chaff  is  set  to  40dB,  and  the  average  ship  RCS  is  set  to  37dB.  Figure  23  shows  the 
probability  of  soft-kill  for  this  stressing  scenario.  It  is  seen  that  despite  the  short  range  where  the 
missile  appears,  chaff  solutions  would  still  protect  the  ship  relatively  well  for  most  wind 
directions. 


Missile  • 
1 


Speed  300  m/s 
range  4,700  m 


Wind  direction 


Figure  23:  Trajectories  of  all  solutions  that  exceed  minimum  probability  of  soft-kill 
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4  Future  Work 


The  SISWS  chaff  solution  algorithm  was  improved  as  described,  in  this  document,  by  adding  a 
monopulse  antenna  model  to  the  seeker,  and  by  using  a  more  physical  approach  to  the  selection  of 
solutions.  The  generated  chaff  solutions  should  provide  more  accurate  responses  for  optimal 
survivability. 

The  next  step  would  consist  of  analysing  the  efficiency  of  the  chaff  solutions  in  SADM  and  also 
comparing  the  results  with  the  SISWS  chaff  solution  algorithm.  After  that,  the  algorithm  can  be 
modified  further  to  include  the  tactic  of  combining  Continuous  Wave  (CW)  jamming  with  chaff 
solutions. 
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Annex  A  Distance  traveled  in  function  of  wind  speed 
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Figure  24:  Arc  of  a  4°  beamwidth  at  range  R. 
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Annex  B  Monopulse  tracking  for  two  RCS  targets 


In  this  annex,  the  monopulse  response  (A/£  channel)  is  analyzed  to  determine  how  well  a  target 
moving  away  from  another  stationary  target  will  be  tracked.  Figure  25  though  Figure  27  shows 
the  delta  sum  ratio  of  the  monopulse  tracker  when  two  targets  of  same  RCS  values  are  separating 
in  time.  Target  1  is  fixed  and  is  located  at  0°.  Target  2  is  moving  away  from  target  1  in  discrete 
steps  from  one  figure  to  the  next.  Red  dots  represent  the  position  of  the  two  targets  at  each  time 
snapshot.  On  Figure  25a,  the  monopulse  response  will  cause  the  tracker  to  be  centered  on  the 
equivalent  RCS  center.  As  target  2  moves  away  from  target  1,  the  monopulse  response  will  cause 
the  tracker  to  follow  the  equivalent  RCS  center,  as  shown  on  Figure  25b.  Around  the  last  track 
position,  the  monopulse  response  will  become  flat,  as  shown  on  Figure  26a.  As  a  result,  the 
monopulse  tracker  will  stay  at  the  last  track  position.  As  target  2  continues  to  move  away,  the 
monopulse  response  will  revert  in  polarity.  This  means  that  the  tracker  will  now  track  the  fix 
target,  as  shown  on  Figure  26b  and  Figure  27,  causing  the  monopulse  to  track  the  fix  target. 
Looking  at  other  values  of  RCS  for  target  2,  it  was  found  that  the  equivalent  RCS  center  will  be 
tracked  until  the  monopulse  system  switches  to  one  of  the  targets.  The  monopulse  tracker  will 
come  back  to  the  ship  whenever  the  ship  RCS  is  greater  or  equal  to  the  chaff  RCS.  Also,  the 
probability  of  recapture  by  the  ship  decreases  as  the  chaff  RCS  exceeds  the  RCS  of  the  ship. 
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Figure  25:  Monopulse  delta  sum  ratio  for  2  targets  as  a  function  of  monopulse  direction 
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Figure  26:  Monopulse  delta  sum  ratio  for  2  targets  as  a  function  of  monopulse  direction 
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Figure  27:  Monopulse  delta  sum  ratio  for  2  targets  as  a  function  of  monopulse  direction 
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ASM 

Antiship  Missile 

CW 

Continuous  Wave 

DND 

Department  of  National  Defence 

DRDC 

Defence  Research  &  Development  Canada 

DRDKIM 

Director  Research  and  Development  Knowledge  and  Information 
Management 

RCS 

Radar  Cross  Section 

RGW 

Range  Gate  Width 

RTG 

Range  To  Go 

R&D 

Research  &  Development 

SADM 

Ship  Air  Defence  Missile 

SISWS 

Shipboard  Integration  of  Sensor  and  Weapons  Systems 

SNR 

Signal  to  Noise  Ratio  (SNR) 

TDP 

Technology  Demonstration  Project 

TTG 

Time  To  Go 
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